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ABSTRACT 

Galaxy formation models typically assume that the size and rotation speed of galaxy disks 
are largely dictated by the mass, concentration, and spin of their surrounding dark matter 
haloes. Equally important, however, are the fraction of baryons in the halo that collect into 
the central galaxy, as well as the net angular momentum that they are able to retain dur- 
ing its assembly process. We explore the latter using a set of four large cosmological N- 
body/gasdynamical simulations drawn from the OWLS (Overwhelmingly Large Simulations) 
project. These runs differ only in their implementation of feedback from supernovae. We find 
that, when expressed as fractions of their virial values, galaxy mass and net angular mo- 
mentum are tightly correlated. Galaxy mass fractions, nrid = M g ai/M V ir, depend strongly 
on feedback, but only weakly on halo mass or spin over the halo mass range explored here 
(M v ; r > 10 11 /i _1 M ). The angular momentum of a galaxy, expressed in units of that of 
its surrounding halo, jd = Jgal/ ■Adr> correlates with rrid in a manner that is insensitive to 
feedback and that deviates strongly from the simple jd = rrid assumption often adopted in 
semi-analytic models of galaxy formation. The rrid-jd correlation implies that, in a given halo, 
galaxy disk size is maximal when the central galaxy makes up a substantial fraction (~ 20- 
30%) of all baryons within the virial radius (i.e., rrid ~ 0.03-0.05). At z = 2, such systems 
may host gaseous disks with radial scale lengths as large as those reported for star-forming 
disks by the SINS survey, even in moderately massive haloes of average spin. Extended disks 
at z = 2 may thus signal the presence of systems where galaxy formation has been particularly 
efficient, rather than the existence of haloes with unusually high spin parameter. 

1 INTRODUCTION One important and robust corollary of these results is that 

the characteristic density of haloes should broadly track the mean 
Galaxy disks are widely assumed to form as gas cools and flows hy- , .. c ., , c ■,, .., , 

J J b 3 density ot the universe and, therefore, increase rapidly with red- 

drodynamically to the center of dark matter haloes. There, baryons , c . r\ a- i a i i in 

J J • j shift. On dimensional grounds, since circular velocity scales like 

settle into thin, rotationally-supported structures whose size and , A rFTTTTn rFt — rv> a v 

' ■> V on WGA1/H oc -v/Gpit , an increase in density, p, implies 

rotation speed are determined by the detailed mass profile of the , ■ .. rnr, ■.. T . r- - r ., , 

v also an increase in the ratio v/H with z. Therefore, it the size and 
system and the net angular momentum of the cooled baryons. In . .. , r , , . , a . ,. , ■ 

f 6 I — s i -i rotation speed ot galaxy disks simply reflect the radius and circu- 

this formation scenario, first laid out by Fall & Efstathiou (1980) 



' i 1 I I ^ ' ^ lar velocity of their surrounding haloes, the size of disks of given 

and worked out in detail by Mo et al. 1 1998. hereafter MMW98). t t - tAt , . ... . 

3 ' 1 1 ^ rotation speed is expected to decrease steadily with increasing z. 

the tight correlations observed to link the mass, size, and rotation 

speed of galaxy disks are thought to reflect analogous correlations The recent discovery of a population of extended disk galaxies 

between the mass, size and spi n of their surrounding dark haloes at 2 = 2 ty the S1 ^ S survey, has confounded this simple expec- 

( see e g iNavarro & SteinmetzkoOOl) tation dGenzel et alj|2006t iFQrster Schreiber et alj|2006l) . In some 

' For dark haloes, these scaling laws have been extensively stud- cases ' these § alaxies have rotation s P eeds comparable to those of 

ied through cosmological N-body simulations, and are now fairly L * dlsks toda y and are as extendedas their z = counterparts, 

well understood. Mass and radius scale so that galaxy haloes are all a result that has led| Bouche et alj 120071) to argue that haloes with 

systems of roughly similar overdensity that, when properly scaled, unusually large spin-much greater than expected in the prevail- 



follow a mass profile well approxim ated by a simple "universal" in 8 ACDM P a r a digm-may be required to explain disks as large as 

fitting formula dNavarro et aljfl99l . [l997l) . Similarly, the net spin observed in the SINS survey. 

of dark haloes, when expressed in dimensionless form (A), appears We note, however, that disk size depends not only on the prop- 

on average to b e independent of mass and/or redshift (see, e.g., erties of dark haloes, but also on the fraction of baryons that col- 

iBett et alj|2007t) . lapse to form the disk, and on the net angular momentum that these 
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Figure 1. Top: Galaxy mass (m^ = M ga x/M v i T ) versus angular momen- 
tum (jd = Jg a i/J v i r ) of simulated galaxies at z = 2. Dots correspond to 
individual systems in the WF2 run; the thick solid curve traces the me- 
dian jd as a function of m^. Other curves show the median for the 
other three runs. The jd- m d relation departs significantly from the often 
assumed jd = relation (top thin solid line). Note that the jd- m d rela- 
tion is approximately the same in all runs and is well approximated by Eq. 
(1): j d = 9.71 m d (l - cxp[-l/(9.71m d )]) (dotted line). The shaded re- 
gion on the right indicates where exceeds the universal baryon fraction, 
/bar ■ Bottom: Distribution ofrrij for each run. We show the histograms for 
WF4 and WF2Dec, our two extreme feedback models, as well as fits with 
Eq.[5]for all four runs. 



lations use a substantially modified version of the s moothed particle 
hydrodynamics code GADGET- 3 dSpringell20oH) . 

We select for this study four OWLS runs that differ from each 
other only in the way energetic feedback from massive stars is im- 
plemented. All these runs assume a standard ACDM cosmology 
consistent with the WMAP-3 results (fi M = 0.238, Ocdm = 
0.1962, n A = 0.762, n b ar = 0.0418, h = 0.73, <t 8 = 0.74, 
n = 0.951) and follow the evolution of 512 3 dark matter and 512 3 
gas particles in a 25 hT 1 Mpc box to z = 2. The mass per baryonic 
particle is ~ 1 x 10 6 Ii~ 1 Mq and 4.7 times higher for the dark mat- 
ter component. The runs adopt a gravitational softening that never 
exceeds 0.5 /i -1 kpc (physical). 

Gas is allowed to cool radiatively using the element-by- 
element implementation of Wiersma et al. (2009a). Gas with den- 
sities exceeding nu = 0.1 cm" 3 is allowed to form stars at 
a pressu re-dependent rate consistent w ith the Schmidt-Kennicutt 
relation JSchave & Dalla Vecchial l2008). The timed release of 1 1 
chemical elements by evolving stellar populations is implemented 
as described in Wiersma et al. (2009b). Feedback energy from su- 
pernovae is incorporated assuming that a fixed fraction of the avail- 
able energy (40% of the total energy released by supernovae for the 
adopted Chabrier IMF) is channeled into "winds" outflowing from 
regions of active star formation. All four runs assume that the to- 
tal energy (per solar mass of stars formed) invested in the outflow 
is the same, but differ in their numerical implementation, which is 
controlled by two parameters: the wind velocity (« w ) and the mass 
loading (77) factor. The parameter 77 specifies the number of gas par- 
ticles among which the feedback energy from a single star particle 
is split, wherea s v w characterizes the outflow velocity of particles 
in the wind (see lDalla Vecchia & Sch ave 2008 for details). 



baryons retain during the disk assembly process. These parameters 
are difficult to estimate, as they likely depend in a complex manner 
on the intricate process of mass accretion, star formation, and the 
regulating effects of energetic feedback that accompany the forma- 
tion of a galaxy. A further complication stems from the fact that 
SINS observations measure the size of the star-forming gaseous 
disk through its redshifted H a emission. This may differ from the 
true spatial extension of all the baryons in the disk, which is what 
theoretical models predict most accurately. 

We address these issues here using a set of large cosmological 
N-body/gasdynamical simulations from the OWLS (Overwhelm- 
ingly Large Simulations) project (Schaye et al, in preparation). In 
this Letter we study the dependence of galaxy disk size on super- 
nova feedback and its consequence for the interpretation of 2 = 2 
galaxy sizes reported by the SINS collaboration. A more general 
analysis of the relation between galaxy mass, angular momentum, 
feedback, and accretion history will be presented in a companion 
paper (Sales et al, in preparation). 



2 THE NUMERICAL SIMULATIONS 



The OWLS project is a collection of over 60 N- 
body/gasdynamical cosmological simulations with varying numer- 
ical resolution and various choices for the implementation of star 
formation, cooling, feedback, and other sub-grid physics. All simu- 



For massive galaxies, the overall effect of feedback, as mea- 
sured by how effectively it regulates star formation and/or removes 
gas from star-forming galaxies, increases with v„ for fixed r\ v J . 
We refer to each of the 4 runs, in order of increasing feedback ef- 
ficiency, as: WF4 (?? = 4 and w w = 424 km/s), WF2 (77 = 2 and 
ti w = 600 km/s), WF1 (77 = 1 and ?j w = 848 km/s) and WF2Dec. 
The latter is equivalent to WF2 but "wind" particles are temporar- 
ily decoupled from the hydrodynamical equations, allowing them 
to freely escape the interstellar medium. 

We focus on all simulated galaxies at the centers of haloes in 
the virial mass range 10 11 < M vir /(/i _1 M ) < 3 x 10 12 . We 
define the virial radius, r v j r , of a system as that of a sphere enclos- 
ing a mean overdensity 178 times the critical value at that redshift. 
All "virial" quantities are measured at or within that radius. Our 
galaxy sample contains ~ 170 objects with between 50, 000 and 
500, 000 particles within the virial radius. The center is given by 
the position of the most-bound particle and the center of mass ve- 
locity corresponds to the average velocity of the 1000 innermost 
particles. Central "galaxies" are identified in each halo with the 
baryonic component of the system defined to be contained within 
r gai = 0.1 r v i r from the halo center. We have explicitly checked 
that this definition includes all of the stars and cold, star-forming 
gas obviously associated with the central object. Since the four 
simulations have identical initial conditions, any difference in the 
galaxy population can be traced directly to the different feedback 
implementation in each run. 
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3 RESULTS 

3.1 Galaxy mass and angular momentum 

The baryonic mass of central galaxies depends strongly on the de- 
tails of the feedback implementation. This is shown in the bottom 
panel of Fig. [JJ where we plot the distribution of galaxy masses 
(gas and stars within r ga i), expressed in units of the virial mass of 
the surrounding halo, = M ga i/M v i r , for our four OWLS runs 
at z — 2. When feedback effects are weak, such as in WF4, most 
of the baryons flow unimpeded to the center of the halo. As a re- 
sult, rrid is on average quite large and approaches in some systems 
the theoretical maximum posed by the universal baryon fraction, 
/bar = fibar/^m = 0.17. At the other extreme, when feedback is 
strong, as in WF2Dec, a large fraction of baryons are pushed out 
of the central galaxy by feedback-driven winds. This leads to low 
values of mj (< 0.02), implying that fewer than ~ 8% of all avail- 
able baryons are retained by the central galaxy in a typical halo. 
Intermediate feedback choices (as in WF1 and WF2) lead to ma 
distributions that straddle those two extremes. In all cases, how- 
ever, the nid distribution is quite broad and insensitive to halo spin 
and mass, at least over the range of masses included in our sample. 

One surprise, given the sensitive dependence of md on feed- 
back discussed above, is the fact that the net angular momentum 
of the central galaxy, when expressed in units of the virial value 
(jd = Jgai/Jvir), correlates with ma, in a manner that is approx- 
imately independent of feedback. This is shown in the top panel 
of Fig.[JJ where the thick line track the median jd as a function of 
md- Points show individual values only for WF2, others are omitted 
for clarity. (The same relation is found in a simulation identical to 
WF2 but with 8 x lower mass resolution, indicating that this result 
is not compromised by numerical resolution.) All four runs follow 
the same md-jd trend, which is well approximated by the relation 



WF2 



3d 



= 9.71 m 2 d (1 - exp[-l/(9.71 m d )}) 



(1) 



where the last factor ensures that jd asymptotically approaches mj 
for large md values (see dotted line in Figure [TJ. Our simulations 
only sample the region 0.005 < mj < 0.1. Care should therefore 
be taken when extrapolating eq. ([TJ beyond these limits. The scatter 
about this relation is not negligible, with an rms of order cr\ og j d ~ 
0.33 about the median trend. Note that this relation implies that 
the specific angular momentum of a galaxy, j m \ — J ga i/Af ga i, 
is quite different from the specific angular momentum of its sur- 
rounding halo, jvir = Jvir/M v i r , unlike what is often assumed in 
semi-analytic models of galaxy formation (J ga i = j v ir is equiv- 
alent to setting jd = md). Indeed, j ml and j vlr are very poorly 
correlated in our simulations ( see also Ivan den Bosch et al]|2002l ; 
button & van den Bosci]|2009t) . 



3.2 Galaxy sizes 



According to the MMW98 formalism, once jd and md are 
specified, the size of a galactic disk may be predicted in a halo 
of known virial radius r v i r , concentration c, and spin parameter 
A. The model assumes that baryons settle onto a thin, rotationally 
supported disk with an exponential radial density profile. The scale 
radius Rd of such a disk is given by 

1 \ja 



R. 



■ fn(jd,m d ,c, A)r vir 



(2) 



(2/ c ) 1/2 rn d 

where /u and f c are known functions that account for the effect of 
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Figure 2. Half mass radii predicted by the MMW98 formalism (a>axis) 
versus the half-mass radii of the stellar (starred symbols) and gaseous (open 
squares) components of simulated galaxies at 2 2. Note the roughly 
constant offset in size between gas and MMW98 predictions. 



halo concentration as well as its adiabatic "contraction" in response 
to the assembly of the galaxy. 

It is instructive to see how this simple prediction fares in com- 
parison with simulated galaxies. This is shown in Fig. [2] where 
we compare half-mass radii measured for galaxies identified in the 
WF2 run with those computed using Eq. l(2j. We use half-mass radii 
instead of exponential scale lengths because it is a more robust size 
estimate, especially when profiles deviate from simple exponen- 
tials, typically as a result of the presence of a spheroid. Starred 
symbols correspond to the stellar component and open squares to 
the gaseous (star forming) disk component. Half-mass radii are 
computed in projection, after using the angular momentum of the 
galaxy to turn the system "face-on". The shaded area indicates radii 
smaller than the gravitational softening. 

On scales exceeding the softening radius, the agreement be- 
tween the size of the stellar component and the MMW98 predic- 
tion is reasonable, especially considering the many simplifying as- 
sumptions in the model. For example, MMW98 assume that all of 
the mass is in an exponential disk, whereas many of the simulated 
galaxies show prominent spheroidal components. 

Encouragingly, the gaseous disks show a much better corre- 
lation with MMW98 model, and seem to differ from R™ MWSS by 
an approximately constant factor n r ~ 1.8. This implies that the 
predictions of Eq. ® may be simply corrected by n r to provide 
reasonably accurate estimates of the size of gaseous disks such as 
those observed in the SINS survey. 



3.3 Disk size and rotation speed 

The MMW98 formalism predicts not only the size of a disk but 
also its rotation speed, V IO %, once the halo parameters (r v i r , c, and 
A) as well as the parameters characterizing the galaxy {md and jd) 
are specified. As discussed in the previous section, although md 
and jd are tightly linked, the mj distribution is very broad, so it 
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is important to explore how much disk sizes may vary solely as a 
result of variations in md- 

We illustrate this in Fig.[5]for the case of a halo with virial ve- 
locity Vvir = 1 50 km/s, c = 4 (the average conce ntration for that 
mass at z = 2; lGao et al.lboOSl ; iDuffv et afj|2008h . and A = 0.04 
(roughly the m edian spin parameter of ACDM haloes; see, e.g., 
iBett et al.ir2007t) . We use the MMW98 formalism to compute the 
size and rotation speed of a gaseous disk galaxy in this halo as a 
function of m,d, assuming that jd(md) follows the median relation 
for run WF2 (thick solid line in Fig.[]} and including the K r correc- 
tion. The result is shown by the dashed curve in Fig.[3] which spans 
a very wide range, from md = 0.01 at the bottom left endpoint to 
md = 0.1 at the upper right. 

A few points are highlighted by this exercise. The first is the 
very wide range in disk size expected from a broad md distribution: 
the disk size changes by a factor of roughly ~ 7, from ~ 600 pc to 
almost ~ 4 kpc (physical), when m<j increases from 0.01 to 0.05. 
This implies that the size of galactic disks is as sensitive to md as 
to the spin parameter. 

The second point is that disk size is maximized for md ~ 
0.05, and declines when mj increases further. This is because 
when md exceeds a certain value, the disk becomes so massive 
that it begins to significantly steepen the central potential, which 
in turn implies that higher rotation speeds are required to achieve 
centrifugal equilibrium. As a result, disks rotate faster but actu- 
ally become smaller, leading the dashed curve in Fig. [3] to decline 
for md > 0.05. Accurate predictions of the size of galactic disks 
clearly require that md and jd be very tightly constrained. 

The pentagons, tria ngles, and circ l es in F i g.[3]show SINS dat a 
fro m the compilations of | Bouche et all ( 120071) . Cresci et al. (2009), 
and Forster Schreiber et al. 1 20091) . Galaxy radii from lBouche et al.l 
torn are H a half-light radii (taken from their Fig. 1) whereas for 
Cresci et al (2009) the radii plotted are the "disk scale lengths" 
derived from their dynamical modeling (taken from their Table 2). 
Radii for the Forster Schreiber et al galaxies are also H a half-light 
radii, as given in their Fig. 18. 

Although we show predicted half-mass gas radii in Fig. [3] we 
have also estimated in our simulated galaxies H a half-light radii by 
weighting each gas particle by its star formation rate and by apply- 
ing an extinction correction based on the metal column through the 
star-forming gas. This procedure gives results which are similar to 
the simple half-mass radii of the gaseous disks measured directly 
from the simulations, suggesting that the comparison of SINS disk 
sizes with the half-mass gas radii of simulated galaxies is fair. 

The dashed-curve in Fig.[3]thus makes clear that even a mod- 
erately massive halo of average concentration and spin may host a 
disk as large as some of those reported in the SINS survey (~ 3-4 
kpc if md ~ 0.03-0.05). Of course, larger disks may also exist; 
for example, in more massive haloes, or in haloes with larger-than- 
average spin, or in systems that scatter above the mean jd-md re- 
lation shown in Fig.[T] For example, the curves in Fig. [3] labelled 
A = 0.04 and A = 0.06 indicate where disks would lie if md is 
fixed at 0.05 and the halo mass is varied. Disks as large as ~ 10 
kpc are clearly allowed for the right combination of A and mj. 

We conclude that the presence of some extended gaseous disks 
at z = 2 is not necessar ily in conflict with the ACDM scenario (see 
also iDekel et alj 2009). However, without strong and tight con- 
straints on the md distribution at z — 2 it is quite difficult to assess 
whether the existence of extended disks in the SINS survey presents 
a real challenge to the ACDM paradigm. Constraining the dis- 
tribution is, unfortunately, a non trivial task where direct simula- 




100 200 
V rol [km/s] 

Figure 3. The radius- velocity plane of observed high redshift galaxies from 
the SINS survey: Bouche et al. (2007, magenta pentagons), Cresci et al. 
(2009, green triangles) and Forster Schreiber et al. (2009, dark circles) (typ- 
ical errors are indicated on the top left). These are compared with the half- 
mass radii of gaseous disks, computed using Eq. {2) and corrected by the 
K r = 1.8 factor. The dashed line corresponds to a disk in a halo with 
Vvi r = 150 km/s, spin parameter A = 0.04, and concentration c = 4; 
with a jd — rtid relation as given by the thick solid curve in Fig [T] Pairs 
of (rrid,jd) values are quoted to help interpretation. The solid curves out- 
line the size- velocity relation expected for disks in haloes with two different 
spin parameters (A = 0.04 and 0.06) and = 0.05. 



tions offer little guidance, given the strong sensitivity of mj to the 
uncertain implementation of feedback highlighted above. 

On the other hand, semi-analytic models suggest that a fairly 
low galaxy formation efficiency is required in order to account for 
the galaxy stellar mass fu nction and its evolution with redshift. 
IConrov & Wechsletl ((2009), for example, argue that fewer than 
~ 10% of baryons should be transformed into stars in galaxy-sized 
haloes. This implies < 0.017. Should this constraint hold true, 
it would imply that very few systems may, on average, reach the 
md ~ 0.05 needed to maximize their disk size. Figure 1 shows that 
many of our galaxies are more massive relative to their haloes than 
inferred by Conroy & Wechsler (2009), particularly if the SN feed- 
back is weak. Unless star formation is quenched in these objects, 
these simulations would therefore overestimate the abundance of 
massive galaxies if they were continued to lower redshifts. 



We can make some progress by assuming that the mj dis- 
tribution, dn/dmd, is independent of halo mass and spin (which 
approximately holds over the mass range considered here) and by 
parametrizing it with a simple function, 



dn/dx — A— j [e 



-((x-i)/o 2 



1]. 



(3) 



where x = mj/Jbar and A is a normalization factor. We note, 
however, that md does depend on mass when extending this analy- 
sis to lower mass halos. Eq. l[3} provides an adequate fit to the md 
distribution in all our simulations, as shown by the thin lines in the 
bottom panel of Fig.Q] The single parameter a x fully specifies the 
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Figure 4. Number density of extended gaseous disks expected in the 
ACDM scenario at z = 2. The different curves show, as a function of the 
median rrid (of a distribution given by Eq. (5)), the number of gaseous disks 
with V ro t > 100 km/s and half-mass radii exceeding 2, 3, and 5 h^kpc. 



distribution; the location of the peak varies monotonically with a x , 
shifting from m d w 0.11 (a x = 0.76) for WF4 to m d » 0.015 
(a x = 0.01) forWF2Dec. 

Adopting Eq. <[3j to describe the ma distribution, and Eq. (T) 
for the jd{md) dependence (including scatter, <J\ gj d = 0.33), 
we can compute the abundance of extended disks expected in 
the ACDM scenario by populating the dark haloes in the simu- 
lations with disks of size consistent with these distributions. For 
each halo we assign, at random, a con centration taken from a log- 
normal distribution centered on c = 4 jGao et alj2008h with width 
(Tio gc = 0.1 dNeto et alj|2007h . and a spin parameter drawn from 
the actual A distribution of the sample. The latter has median 0.028 
and dispersion <7\ og \ = 0.25, roughly consistent with Bett et al 
(2007). 

Fig. [4] shows, as a function of m™ ed //ba, r (median of the 
dn/dx distribution), the number density of gaseous disks with 
Kot > 100 km/s and half-mass radii exceeding 2, 3, and 5 h^ 
kpc (physical) at z = 2. For m™ cd ~ 0.1 /bar, which corresponds 
to the case where ~ 10% of baryons in a galaxy-sized halo have 
been assembled into galaxies, we expect of order ~ 1 extended 
disk per 10 4 cubic comoving Mpc. Even for a m™ ed //b ar as 
low as 0.05, one would still expect of order one such system per 
5 x 10 J cubic h? Mpc. Comparing these numbers to the abun- 
dance of disks at z = 2 from the SINS survey would help to con- 
firm whether the presence of extended disks is indeed in conflict 
with ACDM. 



4 CONCLUSIONS 

We have analyzed the baryonic mass and angular momentum con- 
tent of galaxies identified in four different runs selected from the 
OWLS set of cosmological N-body/gasdynamical simulations. We 
focus on galaxies forming at the centers of haloes with virial mass 



between 10 11 and 3 x 10 12 h~ 1 M Q at z = 2. The four runs ana- 
lyzed differ only in the way the supernova feedback is implemented 
numerically in the simulations. Our main conclusions may be sum- 
marized as follows. 

• The efficiency of feedback, as measured by the fraction of 
baryons that assemble into a central galaxy in each halo, varies 
strongly in our four simulations. In the least efficient feedback run, 
WF4, typically about 50% of the baryons in a halo get assembled 
into the central galaxy. At the other extreme, in run WF2Dec, where 
feedback efficiency is highest, central galaxies are able to retain less 
that 10% of all available baryons. 

• Despite these differences in feedback efficiency, galaxies iden- 
tified in all four simulations follow the same tight relation between 
mass (expressed in units of the halo virial mass, ma) and angular 
momentum (also expressed as a fraction of the virial value, jd)- 
This relation deviates strongly from the simple j d =md assumption 
often used in semi-analytic models of galaxy formation and is well 
approximated by jd = 9.71 (1 — exp[— l/(9.71irid)]), with 
rms dispersion (Tiogj d ~ 0.33. 

• The galaxy mass fraction parameter, rrid, shows a broad dis- 
tribution, with a mean and dispersion that are insensitive to halo 
spin and mass (over the range spanned by our sample). Together 
with the tight jd-rrid relation, this implies that the size of galac- 
tic disks will vary widely, even for galaxies formed in haloes of 
similar mass, concentration, and spin. The size of galactic disks is 
maximized when mj w 0.03-0.05. Such disks can easily match the 
observed size of gaseous disks in the SINS survey, even in moder- 
ately massive haloes of average spin. 

• The abundance of extended disk systems expected at z — 2 for 
ACDM depends critically on the ma distribution. Assuming that 
the form of this distribution may be approximated by Eq. {3), we 
have computed the expected abundance of gaseous disks with half- 
mass radii exceeding Rh = 2, 3, and 5 h^ kpc (physical). If the 
average m d is such that ~ 10% of baryons collect into galaxies, we 
would expect of order ~ 1 system with Rh > 3 kpc in a 10 4 Mpc 3 
(comoving) volume. This predicted abundance could be compared 
with observational estimates, once they become available. 

We conclude that galaxy formation efficiency, through the 
tight jd-rrid relation found in our simulations, plays a crucial role 
in setting the size and rotation speed of galaxy disks. This should 
be taken into account carefully in galaxy formation models. Ex- 
tended disks are expected to exist at z — 2 in the ACDM scenario, 
but in small numbers. The large disks found by the SINS survey 
should not be "typical" galaxies at z = 2, but rather some of the 
largest. The selection of SINS galaxies favors large disks in order 
to enhance the probability of obtaining resolved velocity maps, so 
cannot be considered "typical' either. It seems likely that a num- 
ber of large disks at z = 2 could be accommodated in the ACDM 
scenario without resorting to unexpectedly high halo spins or other 
unusual formation mechanism. 
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